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ELASTIC  DISTORTION  OF  THE  HIGH  PRESSURE  COMPRESSIBILITY  BOMBS  AT  30°  C 

by 

1/ 

Ted  C.  Briggs— 

ABSTRACT 

This  report  contains  equations  representing  the  change  of  com¬ 
pressibility  bomb  volume  as  a  function  of  internal  and  jacket  pres¬ 
sures.  Methods  to  correct  for  elastic  distortion  of  the  bombs  are 

included . 

INTRODUCTION 

2/ 

The  original  Burnett  compressibility  method  (4)~  depended  upon 

V1  +  v2 

the  supposition  that  the  ratio  — r—  of  the  isothermal  volumes  of 

1 

two  containers  is  constant.  In  practice,  the  isothermal  volume  of  a 
container  is  a  function  of  the  confined  pressure,  and  of  the  pressure 
surrounding  the  container. 

The  Thermodynamics  Section  of  the  Helium  Research  Center  is 
using  the  Burnett  method  to  measure  compressibility  factors  of  gases. 
The  compressibility  bombs,  currently  in  use,  are  surrounded  by  oil 

1/  Research  Chemist,  Helium  Research  Center,  Bureau  of  Mines, 

Amarillo,  Texas  _  ,  f 

2/  Underlined  numbers  in  parentheses  refer  to  items  in  the  list  of 

references  at  the  end  of  this  report. 
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jackets.  Pressure  in  the  oil  jackets  can  be  varied  at  will  by  use 
of  an  oil  displacement  pump.  Volumes  of  the  two  compressibility 
bombs  are  functions  of  the  internal  bomb  pressure  and  pressure  in 
the  oil  jackets.  Internal  gas  pressures  can  be  measured  with  a 
high  degree  of  precision  by  use  of  Ruska  Instrument  Corporation 
piston  gage  No.  9274  (72)  .  Jacket  pressures  are  measured  with  a  much 
lower  degree  of  precision  by  use  of  a  Heise  Bourdon  tube  gage.  Pre¬ 
cision  of  measurement  of  the  internal  gas  pressure  ^  0.01  psi  over 
the  pressure  range  30-10,000  psig.  Precision  of  measurement  of  the 
jacket  pressure  is  about  +  5  psi  over  the  pressure  range  0-5,000  psi. 

Elastic  distortion  of  the  PVT  bombs  results  in  an  error  in  the 
experimental  compressibility  factors  unless  a  correction  is  made  for 
distortion  effects.  The  purpose  of  this  report  is  to  develop  equa¬ 
tions  to  represent  the  change  of  bomb  volumes  V^  and  V^  +  V^  as  a 
function  of  the  internal  and  jacket  pressures,  and  to  present  methods 
to  correct  for  elastic  distortion  of  the  bombs. 

DESIGNATION  OF  THE  COMPONENT  VOLUMES  OF  V^^  AND  V2 
A  Burnett  apparatus  consists  of  two  containers  designated  as  V^ 
and  V2 .  The  original  method  consisted  of  a  series  of  expansions  of 
gas  from  V.  into  an  evacuated  V„ .  To  develop  equations  to  represent 
V  and  V-  +  V„  as  a  function  of  the  internal  and  jacket  pressures, 

JL  JL 

one  must  know  the  dimensions  of  V^  and  V2 ,  and  the  magnitude  of  the 
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jacketed  and  unjacketed  portions  of  and  For  the  following 

derivations  notations  are  assigned  to  the  various  volumes. 

V?  =  volume  of  V,  at  zero  internal  and  jacket  pressures 

11 

V°  =  volume  of  at  zero  internal  and  jacket  pressures 

V°  =  volume  of  jacketed  bomb  portion  of  at  zero 
.1  internal  and  jacket  pressures 

V°  =  volume  of  jacketed  bomb  portion  of  at  zero 
2  internal  and  jacket  pressures 

V°  =  volume  of  unjacketed  tubing  portion  of  at 

1  zero  internal  pressure 

V°  =  volume  of  unjacketed  tubing  portion  of  at 

2  zero  internal  pressure 

V°  =  volume  of  unjacketed  fittings  connected  to 

1  at  zero  internal  pressure 

V°  =  volume  of  unjacketed  fittings  connected  to 

2  at  zero  internal  pressure 

CHANGE  OF  JACKETED  BOMB  VOLUME  AS  A  FUNCTION  OF 
INTERNAL  AND  JACKET  PRESSURES 

The  PVT  bombs  were  machined  from  303  stainless  steel  bar  stock. 
Dimensions  of  one  of  the  jacketed  bombs  are  listed  in  figure  1.  In 
the  following  derivations  a  number  of  simplifying  assumptions  are 
made.  Important  assumptions  are  numbered. 

o  o 

Assumption  1:  Dimensions  of  jacketed  bomb  volumes  V,  and  V, 

12 

are  identical. 

From  the  dimensions  of  figure  1,  the  jacketed  bomb  volume  is 

3 

estimated  to  be  4.649  in  .  Examination  of  figure  1  reveals  that  the 
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FIGURE  I. -Jacketed  High-Pressure  Compressibility  Bomb 
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jacketed  bomb  volume  consists  of  several  volumes  of  various  wall 
thickness  and  of  various  shapes.  It  would  be  difficult  to  compute 
the  exact  pressure  function  of  each  of  the  various  jacketed  volumes; 
therefore,  a  simplifying  assumption  is  made. 

Assumption  2:  The  jacketed  volume  consists  of  a  uniform, 
thick-wall,  closed  plane  end  cylinder,  with  an  internal  radius  of 
0.5  inch,  a  wall  thickness  of  0.75  inch,  and  an  internal  volume  of 
4.649  in^ . 

Assumption  3°.  Distortion  of  the  internal  bombs,  tubing,  and 
fittings  is  elastic. 

A  number  of  references  are  available  in  the  literature  per¬ 
taining  to  elastic  distortion  of  cylinders.  Timoshenko  (8)  presents 
an  equation  for  the  change  of  radius  of  an  open-end  cylinder.  Love 
(6)  developed  an  equation  for  the  change  of  volume  of  a  thick-wall 
closed-end  cylinder.  Newitt  (_7)  derived  equations  for  the  change  of 
radius  and  length  for  both  open-end  and  closed-end  thick-wall  cylin¬ 
ders.  These  references  are  used  to  calculate  the  change  of  volume 
of  the  internal  bombs,  tubing,  and  fittings.  Volume  of  a  cylinder 
is  given  by  equation  (1)  . 

V  =  tt  r^L  (i) 

Take  the  derivative  of  equation  (1)  with  respect  to  the  radius  and 
length,  and  replace  the  derivatives  with  delta  quantities  to  obtain 
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equation  (2)  . 

AV 


2tt 


LAr+rrr  AL 


(2) 


Divide  equation  (2)  by  equation  (1)  to  obtain  equation  (3) 


AV  =  2Ar  AL 
V  r  L 


(3) 


A  jacketed  cylinder  with  closed  ends  is  subject  to  a  longi¬ 
tudinal  stress  due  to  pressure  acting  on  the  ends.  The  elongation 
is  given  by  equation  (4)  ,  references  (6}  and  (7)  . 


AL 

L 


( 1-2ct) (a2P .  -  b2P.) 

i _ ]_ 

2  2 
E(bZ  -  aZ) 


(4) 


Change  of  internal  radius  of  a  thick-wall,  closed-end,  jacketed 
cylinder  is  given  by  equation  (5)  ,  reference  (_7)  . 


Ar 

r 


(1-2ct) (a2P .  -  b2P.) 

i _ jl_ 

2  2 
E(bZ  -  aZ) 


(1+ct) (P1  -  P.)b2 

-,.2  2n 

E(b  -  a  ) 


(5) 


Substitute  equations  (4)  and  (5)  into  equation  (3)  to  obtain  equation 

(6). 


AV, 


r 


V 


o 


_ 2  2n 
E(bb  '  V 


3 (1-2ct) (a,  P. 

.  b  i 


b2P.)  +  2(l+a)(P  -P  )b2 

b  J  l  J  b 


n  3/ 


(6) 


3/  The  subscript,  b  denotes  bomb  dimensions. 
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The  terms  in  the  preceding  equations  are: 


V 

L 

AL 


b 


cylinder  volume 

internal  cylinder  length 

change  of  internal  cylinder  length  due 
to  pressure  distortion 

internal  cylinder  radius 

change  of  internal  cylinder  radius 
due  to  pressure  distortion 

external  cylinder  radius 


P  =  internal  cylinder  pressure 
i 

P.  =  jacket  pressure 
CT  =  Poisson's  ratio 


E  =  Young's  modulus 

Equation  (6)  can  be  rearranged  to  give  equation  (7) . 


3 ( 1-2ct) a2  +  2 (1+ct) b2 

E<bb  -  3b> 


P. 

i 


(5-4ct) b2 
E(bb  "  ab>  ^ 


(7) 


Equation  (7)  is  in  the  form  of  equation  (8) . 


c/P.  -  p'p. 
i  J 


(8) 


where 

3 (1-2ct) a2  +  2 (1+a) b2 
a 


E(b2b  -  a2) 


(9) 


/ 


■ 

'  . 
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and 

(5-4cr)tv 

p'  -  -i — r 

E<bb  -  V 

Equation  (10)  can  be  rearranged  to  give  equation  (11) . 


(5-4cr)b3 


(10) 


(ID 


CHANGE  OF  UNJACKETED  TUBING  VOLUME  AS  A 
FUNCTION  OF  THE  INTERNAL  PRESSURE 

The  valves,  fittings,  and  bombs  of  and  V ^  are  connected  by 
1/8  inch  od,  0.035  inch  wall  thickness,  high  pressure,  stainless 
steel  tubing.  The  connecting  tubing  constitutes  a  portion  of  the 
unjacketed  volumes  of  and  Figure  2  is  an  unsealed  block  dia¬ 

gram  showing  the  components  of  and  •  The  total  length  of  tubing 
in  V2  is  20.5  inches,  and  the  total  length  of  tubing  in  is  13.6 
inches . 

V°  =  0.032  in3 

C1 

V°  =  0.049  in3 

t2 

The  change  of  tubing  radius  is  calculated  from  equation  (12) , 
reference  (7) ,  page  45.  Equation  (12)  gives  the  change  of  internal 
radius  of  a  thick-wall,  closed-end,  unjacketed  cylinder  as  a  function 


. 

*  ...  .  I  .  v"  !  8\| 


FIGURE  2. -Block  Diagram  of  Volumes  V(  and  V ^ 
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of  the  internal  pressure. 


Ar 

r 


(l-2a)a2P.  (i+a)b2P. 

1  1 

9  9  +  9  9 

E(b  -  a  )  E(b  -  a  ) 


(12) 


The  change  of  tubing  length  is  calculated  as  if  the  tubing  were 
a  closed-end  cylinder.  The  change  of  tubing  length  as  a  function  of 
internal  pressure  is  calculated  from  equation  (13)  ,  reference  (7) , 
page  45 . 


AL 

L 


(1-2ct) a2P . 
_ _ i 

2  2 

E(b  -  a  ) 


(13) 


Equations  (12)  and  (13)  are  substituted  into  equation  (3)  to  obtain 
equation  (14) . 


+  2(l+CT)b2 


4/ 


(14) 


4/  The  subscript  t  denotes  tubing  dimensions. 


The  change  of  tubing  volume  as  a  function  of  the  internal  pressure 
was  also  calculated  as  if  the  tubing  were  an  open-end  cylinder.  In 
one  case,  the  change  of  tubing  volume  was  calculated  for  an  open-end 
cylinder  with  the  assumption  that  the  change  of  volume  is  due  to  a 
change  of  radius  only;  and  in  another  case,  with  the  assumption  that 


13 


the  volume  change  is  due  to  a  change  of  radius  and  a  decrease  in 

AVt 

length.  The  various  assumptions  give  different  values  for  — —  ; 


A V 


V 


however,  the  effect  of  the  different  values  for  — —  on  the  change 


V 


o 


of  and  with  pressure  is  negligible. 


CHANGE  OF  UN JACKETED  FITTING  VOLUME  AS  A 
FUNCTION  OF  THE  INTERNAL  PRESSURE 


It  would  be  difficult  to  estimate  the  exact  pressure  function 

of  the  fittings  connected  to  and  V  .  The  fittings  which  consti- 

o  3 

tute  V,.  are  the  Ruska  differential  pressure  cell  (0.038  in  ) ,  a 
f 

1 

3  3 

Ruska  tee  (0.014  in  ),  a  Ruska  valve  (0.002  in  ),  and  two  Ruska 

connectors  (0.022  in^) .  The  fittings  which  constitute  V°  are 

2 


3  3 

three  Ruska  valves  (0.006  in  )  ,  three  Ruska  connectors  (0.033  in“ ) , 

3 

and  a  Ruska  cross  (0.012  in  ). 

O 

V°  =  0.076  in 

rl 


V 


0.051 


.  3 
in 


Assumption  4°.  (a)  The  change  of  fitting  volume  can  be  com¬ 

puted  as  if  the  fittings  were  a  thick-wall,  open-end  cylinder,  with 
an  id  of  0.070  inch,  and  an  od  of  0.32  inch.  (b)  The  change  of 
length  of  the  hypothetical  cylinder  is  negligible. 


. 
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The  change  of  fitting  volume  is  calculated  from  equation  (15) . 


AV 


f 


V 


o 

f 


2(1-ct)  +  2(l+cr)b^ 

2  2 

E(bf  “  af ) 


1/ 


P. 

l 


(15) 


5/  The  subscript  f  denotes  hypothetical  fitting  dimensions. 


The  assumed  fitting  dimensions  are  a  rough  estimate  of  the  average 
fitting  dimensions;  therefore,  calculation  of  the  change  of  fitting 
volume  with  pressure  is,  at  best,  approximate. 

EXPERIMENTAL  DISTORTION  COEFFICIENT 
AND  CALCULATION  OF  YOUNG'S  MODULUS 

The  volume  of  as  a  function  of  jacket  and  internal  pressures 

can  be  written  as  equation  (16)  ,  reference  (_3)  . 

V,  =  V?(l  +  aP„  -  PP.)  (16) 

1  1  i  J 

Equation  (16)  can  be  written  as  equation  (17)  when  AV^  =  -  V^. 


AV 

—  =  a P.  -  3P 

V° *  1  j 

1 


(17) 


The  coefficient  3  was  determined  to  be  1.5163  x  10  +  0.0068  x  10  , 

psi"1  at  30°  C  (3).  The  change  of  due  to  change  of  jacket  pres¬ 
sure  is  essentially  equal  to  the  change  of  the  jacketed  bomb  volume 


. 


- 


.  ....  .  •  ' 


*' 
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due  to  change  of  jacket  pressure 


3=3 


v 


v 


.0 


V.  =  4.649  in' 

bi 


V°  =  4.757  in3 


3 '  =  1.5515  x  10”7,  psi"1 


Young's  modulus  can  be  calculated  from  equation  (11) 


(5  -  4a) b 


E  = 


b 


3 7 (bb  ab> 


(ID 


a  =  0.305  (I) 


a,  =  0 . 5  in 

b 


b,  =  1 . 25  in 

b 


E  =  29.004  x  10  psi 


CHANGE  OF  VOLUMES  Vx  AND  V ^  +  V2  AS  A 
FUNCTION  OF  THE  INTERNAL  AND  JACKET  PRESSURES 

Equation  (7)  represents  the  change  of  jacketed  bomb  volume  as 

a  function  of  the  internal  and  jacket  pressures. 

AVb  3(l-2a)ab  +  2(l+CT)bj;  (5-4a)t>b 

Vb  E^bb  -  3b>  E(bb  -  ab>  Pj 
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Numerical  values  for  the  constants  are  substituted  into  equation  (7) 
to  obtain  equation  (18). 

a,  =  0.5  in 

b 

b  =  1.25  in 

b 

ct  =  0.305 

E  =  29.004  x  10^  psi 

=  1.1481  x  10“7  P.  -  1.5515  x  10"7P.  (18) 

i  J 

(14)  expresses  the  change  of  unjacketed  tubing  volume 
of  the  internal  pressure. 


A4 


Equation 
as  a  function 


AV  3(l-2cj)a^  +  2  ( 1+cr)  b^ 

\T°  V/'K2  2\ 

\  E(bt  '  V 


(14) 


Assumption  5:  Young's  modulus  and  Poisson's  ratio  for  the 
bombs,  tubing,  and  fittings  are  identical. 

Numerical  values  for  the  constants  are  substituted  into  equa- 
tion  (14)  to  obtain  equation  (19). 

a  =  0.028  in 

b^  =  0.063  in 

a  =  0.305 


■ 


. 


■ 


17 


E  =  29.004  x  10  psi 


A V. 


V 


o 


1.2206  x  10'7  P. 

l 


(19) 


Equation  (15)  represents  the  change  of  unjacketed  fitting  volume 
as  a  function  of  the  internal  pressure. 


AV 


V 


o 


2(l-a)a2  +  2(l+cr)b2 
E(b2f  -  a2) 


(15) 


The  constants 


a,.  =  0.035  in 


b =  0.160  in 


=  0.305 


£ 

E  =  29.004  x  10  psj. 


are  substituted  into  equation  (15)  to  obtain  equation  (20) 


AV 


V 


o 


f  =  0.9692  x  10"7P 


(20) 


V°  =  V°  +  V°  +  V° 

1  bl  tl  fl 


(21) 


V°  =  V°  +  V°  +  V° 

2  b2  b2  f2 


(22) 


: 


.  f 


r  or  -  ' . 


- 


.. 


••  ■  •  •  j 


•  ■  V 
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A V,  =  AV,  +  AV  +  AV 

1  bi  4  fi 


(23) 


AV„  =  AV,  +  AV  +  AVC 
2  b2  t2  f2 


(24) 


AV. 


V 


o 


AV,  AV  AV 

bi  4  fi 

— i  +  — =•  +  — =■ 

V°  v°  v° 


(25) 


A(V1  +  V 

V°  +  V° 

1  2 


A(V  +  v  )  A(V  +  V  )  A(V  +  V  ) 

bl  b2  4  2  fl  4 

- +  - ± - £_  +  - 1 - —  (26) 


v°  +  v° 

1  2 


v°+v° 


v°  +  v° 

1  2 


Volumes  estimated  from  the  component  dimensions  are: 

V?  =  4.649  in3 


V°  =  4.649  in3 

b  „ 


V°  =  0.032  in3 

4 


V°  =  0.049  in3 

4 


Q 

V°  =  0.076  in 

4 


V°  =  0.051  in3 

4 


V°  =  4.757  in3 


- 


:■ 

■ 


. 


■ 


-  , 


. 
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V°  =  4.749  in3 


V°  +  V° 
1  2 


9.506  in3 


o  ,  „o 


From  equation  (18)  and  the  dimensions  of  and 


AV, 


V 


o 


o 


AVr  V 

h  .  n 


V 


o 


V 


o 


=  1.1220  x  10_7P.  -  1.5163  x  10_7P.  (27) 

i  J 


o  ,  ,,o 


From  equation  (19)  and  the  dimensions  of  and  Vt 


AV. 


V 


o 


AV.  V 


o 


V  V 


o 


0.0082  x  10_7P. 

L 


(28) 


o  ,  „o 


From  equation  (20)  and  the  dimensions  of  and 


AV 


AV 


V 


o 


V 


o 


V 


o 


V. 
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0.0155  x  10"7P. 
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(29) 


Equations  (27)  ,  (28)  ,  and  (29)  are  substituted  into  equation 
(25)  to  obtain  equation  (30) . 


AV 
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-  =  1.1457  x  10_7P.  -  1.5163  x  10_7P. 

o  i  J 


P.  and  P.  are  in  psi 
1  J 


(30) 
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Equation  (30)  represents  the  change  of  volume  V°  as  a  function 
of  internal  and  jacket  pressures. 

Equation  (31)  is  derived  from  equations  (18)  ,  (19)  ,  (20)  ,  (26)  , 
and  the  component  dimensions. 


A(V1  +  V  -7  -7 

— - - —  =  1.1463  x  10  'p.  -  1.5176  x  10  P.  (31) 

V°  +  V°  1  J 

P.  and  P.  are  in  psi 
r  J  r 


Equation  (31)  represents  the  change  of  volume  (V^  +  V°)  as  a  func¬ 
tion  of  internal  and  jacket  pressures. 


Suppose  one  wishes  to  adjust  the  jacket  pressure  so  that  the 

V°  +  V°x 

f  l  2\ 

1°  l  J  will  be  a  constant.  When  the  gas  sample  is 


volume  ratio 


V 


confined  in  ,  the  jacket  pressure  should  be  adjusted  to  0.7556  of 
the  internal  pressure  for  AV^  to  be  zero.  When  the  gas  sample  is 
confined  in  V  +  V2 ,  the  jacket  pressure  should  be  adjusted  to  0.7553 
of  the  internal  pressure  for  A(V^  +  vp  to  be  zero.  Adjusting  the 
jacket  pressure  for  each  measured  internal  pressure  is  one  way  to 
compensate  for  elastic  distortion  of  the  PVT  bombs. 

CORRECTION  OF  THE  VOLUME  RATIO  FOR  ELASTIC  DISTORTION 


Another  method  of  correcting  for  elastic  distortion  of  the  PVT 
bombs  is  the  method  used  by  Canfield  (5) .  Canfield's  method  con¬ 
sists  of  correcting  the  volume  ratio  for  elastic  distortion. 
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,V?  +  V° 


/  1  2\ 

The  volume  ratio  l - - - J  at  zero  pressure  is  defined  to  be 

Vi 

N  .  The  volume  ratio  for  the  r  th  expansion  is  given  by  equation 
(32). 


N 

r 


1  + 


N 


A<vi  +  Vi 

V°  +  V° 

1  2 


at  P 


o 


1  + 


AV 


V 


o 


at  P 


r-1 


(32) 


P^_  ^  is  the  pressure  before  the  r  th  expansion  and  P^  is  the  pressure 
after  the  r  th  expansion. 

r  =  1,2,  3,  4 ,  5... 

Substituting  equations  (30)  and  (31)  into  equation  (32)  gives  equation 
(33). 


(1  +  1.1463  x  iO‘7P.  -  1.5176  x  10_7P.) 

i  j  at  P 

N  =  N  — - - — - - —  (33) 

r  °  (1  +  1.1457  x  10  P.  -  1.5163  x  10  P.)  at  P  , 

l  j  r-1 


Let  the  jacket  pressure  be  zero  and  equation  (33)  reduces  to  equa¬ 
tion  (34)  . 


(1  +  1.1463  x  10~7P.)  at  P 
_ _ i_ _ r_ 

(1  +  1.1457  x  10_7P.)  at  P  , 

i  r-1 


(34) 
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The  fundamental  equation  for  calculation  of  a  compressibility 
factor  by  the  Burnett  method  is  equation  (35) . 


Z 

r 


(35) 


Substitution  of  equation  (34)  into  equation  (35)  gives  equation 
(36). 


Z 

r 


(1  +  k  P  )(i  +  k  P  ) 
No  (1  +  k2Pfl)(l  +  k2P1) 


(I  +  k.P  ) 
_ 1  r 

(i  +  k2pr_l} 


(36) 


where 


k^  =  1.1463  x  10  7,  psi  \  at  30°  C 
k^  =  1.1457  x  10  ^ ,  psi  \  at  30°  C 


A  compressibility  factor,  obtained  by  the  Burnett  method  and  corrected 
for  elastic  distortion  of  the  PVT  bombs,  can  be  calculated  from  equa¬ 
tion  (36)  . 

The  constants  k^  and  k^  are  dependent  upon  the  dimensions  of  V? 
and  ,  Poisson's  ratio  for  and  V and  Young's  modulus  for 
and  V  .  Young's  modulus  is  temperature  dependent;  therefore,  the  con¬ 
stants  k^  and  k^  are  temperature  dependent. 
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